Experimental study of quasi-elastic neutrino interactions on Ar with a 
liquid Ar TPC exposed to the WANF neutrino beam 
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We present results from the first exposure of a liquid Ar time projection cliamber to a neutrino beam. The data 
have been collected in 1997 with a 50 liter ICARUS-like chamber located between the CHORUS and NOMAD 
experiment at the CERN West Area Neutrino Facility. We focus on the analysis of quasi-elastic interactions; 
despite the limited size of the detector, nuclear effects beyond Fermi motion and Pauli blocking have been observed 
as perturbations to the pure quasi-elastic kinematics. 



1. Introduction 

In 1997 a liquid argon time projection cham- 
ber (LArTPC) with a fiducial volume of about 
50 1 was exposed to the multi-GeV wide-band 
neutrino beam of the CERN West Area Neu- 
trino Facility (WANF) P2], during the NO- 
MAD and CHORUS gj data taking. The 
LArTPC was placed on a platform 4.5 m high, 
between the CHORUS and NOMAD detectors. 
The modest size of the LArTPC fiducial volume 
(^^50 liters) made necessary a muon spectrome- 
ter downstream the TPC. A coincidence with the 
NOMAD DAQ was set up to use the detectors lo- 
cated into the NOMAD magnetic dipole as a spec- 
trometer. During the test 1.2 x 10^^ p.o.t. were 
integrated, and about 10^ triggers recorded. The 
test was run jointly by a group from Milano Uni- 
versity and the ICARUS collaboration 0. The 
data accumulated in 1997 still offer the unique 
opportunity of an experimental study on neutrino 
interactions on Ar nuclei. At this time of writ- 
ing, final results from the experiment have not 
yet been reported in any systematic fashion (but 
see |7I8| for preliminary results and a more de- 
tailed description of the experimental apparatus) ; 
a comprehensive set of results from the 1997 test 
is presented here for the first time, while a more 
thorough and final paper is in preparation. 
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2. The experimental setup 

The LArTPC had an active volume of 
32x32x46.8 cm^, enclosed in a stainless steel ves- 
sel in the shape of a bowed-bottom cylinder 90 cm 
high with a radius of 35 cm. The active volume 
contained 67 kg of liquid Ar. Ionization electrons 
produced by the passage of charged particles in 
liquid Ar drift toward the anode under the action 
of 214 V/cm uniform electric field. The readout 
electrodes are two parallel planes of wires run- 
ning orthogonally; the planes are at a distance 
of 4 mm and each plane has 128 wires. Each 
stainless steel wire has a diameter of 100 /um; the 
distance between the wires is 2.54 mm. The ion- 
ization electrons drift through the first (induc- 
tion) wire plane so that the integrated induction 
signal is zero. The collection plane then collects 
the drifting electrons. The mean charge per unit 
pitch (2.54 mm) for a m.i.p. crossing the cham- 
ber parallel to the readout planes, corresponds to 
about 2.3 X 10"* electrons (~4 fC). 

The chamber has been exposed to the v beam 
produced at the CERN WANF. The primary 
450 GeV protons from the CERN SPS were ex- 
tracted every 14.4 s in two spills (1.8 x 10^^ pro- 
tons per spill, on average) of 6 ms duration each 
and separated by 2.5 s, and hit a segmented Be 
target. Secondaries were selected in momentum 
and focused by a system of collimators and mag- 
netic lenses. The experimental area was located 
940 m downstream the target. The mean energy 
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of the reaching the detectors is 24.3 GeV, while 
contaminations from other flavors are below 7% 
for and ~ 1% for Vf. [2]. 

The trigger was provided by a set of plastic 
scintillators, located downstream the chamber be- 
fore the NOMAD apparatus. Incoming charged 
particle were vetoed by 5 large plastic scintillators 
mounted in front of the chamber and by the last 
scintillator plane of CHORUS, which vetoed par- 
ticles deflected by the CHORUS magnetic field 
entering the chamber at large angles with respect 
to the nominal beam direction. 

The local trigger required the coincidence of the 
SPS beam spill, at least one of the trigger scin- 
tillators, and the two trigger scintillator planes of 
NOMAD (Tl and T2 in ^). Moreover, a trigger 
was rejected if the NOMAD acquisition system 
was in BUSY mode or if it came less than 500 /zs 
after the previous trigger, inhibiting the occur- 
rence of overlapping. In this configuration 11% 
of the recorded events had a vertex in the fiducial 
volume of the TPC. The request of a charged par- 
ticle triggering the chamber locally and reaching 
NOMAD up to Tl and T2 hmited the useful data 
to j/^ charged-current interactions. 

3. Event reconstruction and particle iden- 
tification 

Interactions in the TPC fiducial volume were 
fully imaged in two 2D images with a common 
coordinate (time), with full calorimetric informa- 
tion associated with each point. The common co- 
ordinate allows an unambiguous 3D reconstruc- 
tion of the event, at least for simple topologies. 
Each 2D-image represents the signal amplitude 
digitized by the ADC (in a linear gray scale) ver- 
sus the time sample (drift coordinate) and the 
wire number. Fig. ^ shows the image of a quasi- 
elastic (QE) event in the induction and collection 
planes. Units are wire number versus drift time 
in /is. Appropriate filtering techniques have been 
applied to deconvolve the instrumental response 
and precisely associate the charge deposition to 
each pixel in the image. For the collection wires a 
universal response function has been determined 
empirically from events in the m.i.p. sample (el- 
ementary charge deposition). The response was 
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Figure 1. The filtered image of an event recon- 
structed as Vfji [i~Vi with a m.i.p. leaving 
the TPC fiducial volume and an identified stop- 
ping proton. Top induction view; bottom collec- 
tion view. 



then deconvolved by means of a Discrete Fourier 
Analysis to reconstruct the topology of the actual 
charge deposition in the event. A different ap- 
proach was adopted for the induction view, where 
the charge measurement in that view was not ac- 
complished and only the event position, measured 
by the zero-crossing of the induced signal, was re- 
constructed. 

3.1. Proton reconstruction and p/vr sepa- 
ration 

For a proton, the identification and momentum 
measurement were performed using only informa- 
tion from the TPC. The discrimination between 
protons and charged pions is performed exploiting 
the different behavior of energy loss as a function 
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of the range. For the case of candidates stop- 
ping in the fiducial volume of the chamber, as 
the ones of the "golden sub-sample" considered 
in Sec^ the discrimination is highly simplified: 
it is based on the consistency between the recon- 
structed kinetic energy from the range and the in- 
tegrated charge deposited by the candidate along 
its path-length. For partially contained candi- 
dates, a much better assessment of the charge re- 
sponse of the TPC and of quenching effects [TUI 
is required, since proton identification is based 
only on the visible pattern of the energy loss. 
Such events are not considered here. For fully 
contained protons, the momentum uncertainty is 
dominated by the finite pitch of the wires. The 
equivalent pitch in the vertical direction (drift di- 
rection) is much smaller due to the high sampling 
rate of the fast ADC (360 /im). For horizontal 
protons of p = 400 MeV/c the uncertainty is of 
the order of 7 MeV/c. The angular resolution in 
the collection or induction view depends on the 
number of wires N hit by the particle along its 
path; it is cr ~ 0.36\/T2/(2.54 N^/^), correspond- 
ing e.g. to 15 mrad for N=10. 

3.2. Muon reconstruction 

The kinematic reconstruction of the outgoing 
muons exploits the tracking capability of NO- 
MAD. An event triggering the chamber will have 
at least one penetrating track reaching the Tl 
and T2 trigger scintillators bracketing the TRD 
of NOMAD 9 . The corresponding track, nearly 
horizontal at the entrance of the NOMAD drift 
chamber volume, is reconstructed with an aver- 
age momentum precision of (Tp/p ~ 0.05/a/L © 
0.008p/L^/2, L being the visible range in the vol- 
ume itself expressed in meters and p the particle 
momentum in GeV. A 10 GeV horizontal muon 
crossing all the chambers (L ^ 5m) is recon- 
structed with a precision of 2.2%. The recon- 
structed particle is traced back to the TPC ac- 
counting for the magnetic field and the presence 
of the forward NOMAD calorimeter. The latter 
introduces the dominant uncertainty on the muon 
transverse momentum, due to multiple scatter- 
ing (MS) in iron (190 cm for a horizontal muon). 
For small scattering angles {6^1 rad) the MS 
uncertainty on the transverse momentum is inde- 



pendent of p and it turns out to be ~ 140 MeV. 
The correctness of the back-tracing procedure has 
been cross-checked comparing the direction an- 
gles of the particles belonging to the m.i.p. cali- 
bration sample as measured by the TPC with the 
corresponding quantity from NOMAD. 

4. Analysis of quasi elastic i^^ interaction 

4.1. Event selection and rates 

The commissioning of the detector ended in 
August 1997. After that, 1.21 x 10^^ protons 
on target have been integrated. The trigger ef- 
ficiency was monitored during data taking and 
its integrated value is 97%. Additional losses are 
due to the TPC (3%) and NOMAD (15%) dead 
time and to detector faults. The effective live- 
time was 75%. 81,000 events had been recorded 
and analyzed by visual scanning. A minimum- 
bias sample has been obtained requiring no more 
than two tracks exiting the chamber and an ar- 
bitrary number of contained tracks. This sample 
contains the golden sample of QE interactions and 
other samples used to validate the simulation for 
background estimation. The golden sample con- 
sists of events with an identified proton of kinetic 
energy (Tp) larger than 50 MeV fully contained in 
the TPC and one muon whose direction extrapo- 
lated from NOMAD matches the outgoing track 
in the TPC. The distance of the interaction ver- 
tex from any of the TPC walls has to be greater 
than 1 cm. The muon candidate track projected 
onto the wire plane must be longer than 12 wire 
pitches. The event is accepted even in the pres- 
ence of other stopping particles, as far as their ki- 
netic energy (in the proton hypothesis) does not 
exceed the Tp of the leading proton. If tracks 
other than the identified muon leave the TPC or 
at least one converted photon with energy greater 
than 10 MeV is present in the fiducial volume, the 
event is rejected. The tightness of these selections 
defines a very clear topology for visual scanning 
(Fig.UJ. 

The lower Tp cut combined with the request 
of containment is very severe since the chamber 
volume is small. On the other hand, at lower 
Tp the proton range is comparable with the wire 
pitch and neither the proton momentum nor the 
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interaction vertex can be reconstructed with due 
precision. Moreover, for Tp > 50 MeV the 7r*/p 
misidentification probability is negligible. The 
golden sub-sample contains pure QE interactions 
(i/^ n fj,~ p), an intrinsic background domi- 
nated by resonant productions followed by pion 
absorption in the nucleus { v^, p A"*""*"/!^ — > 
/i^ p 7r+ , n — > A+^^ /i^ p i:^) and an 
instrumental background due to unidentified tt'^'s 
{v^ n —> A"*"/!" —> /i~ p tt'^). On the other hand, 
the n —t A+/i^ —t fi^ n 7r+ contamination is 
negligible in this tightly selected sample. 

The efhciency of the selections for QE inter- 
actions and their intrinsic background has been 
evaluated by Monte Carlo experimentation based 
on the FLUKA code ^Jj and it turned out to 
be 17%. Similarly, MC provided the inefhciency 
of the vetoing selections for n A+/i^ 
/i~ p TT*^, i.e. the probability to miss both the de- 
cay photons of the tt*^ —> 77 or the e~^e~j system 
in case of 7r° Dalitz decay. The corresponding 
contamination of the golden sub-sample is esti- 
mated to be 13%, which can be checked by test 
samples extracted from the minimum bias events. 
The test samples consist of golden events with 
one (iVi) or two (A^2) converted photons pointing 
to the interaction vertex. Assuming the gamma 
identification probability to be uncorrelated for 
the two photons, we have A^i = 2N{1 — e^)e-y and 



Ne^, N being the (unknown) overall rate of 
p ii~ 7r° final states and e-y the photon identifica- 
tion efficiency; = (A^i/2A^2 + 1)~^ turned out to 
be (43± 9) %. Hence, the probability of missing 
both gammas is (32±10) % to be compared with 
the MC calculation of 20.4%. Accounting for the 
detector livetime, trigger and selection efficiency 
and the sample purity, we expect 73.5 events and 
we observe 61. 



• Data 

— MC (Fluka) 

— MC (SW) 
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4.2. Analysis of the kinematic distortion 
due to nuclear matter 

In spite of the limited statistics, the golden sub- 
sample provides informations on the basic mech- 
anisms that modify the kinematics of neutrino- 
nucleus interactions with respect to the corre- 
sponding neutrino- nucleon process. Nuclear mat- 
ter perturbs the initial state of the interaction 
through Fermi motion; it also affects the for- 



Figure 2. From top to bottom: distribution of the 
muon energy, muon transverse momentum and 
fi — v angle for the golden sub-sample. The con- 
tinuous red (green) line is the expectation from 
FLUKA (Saxon- Woods) convoluted with the de- 
tector response. 
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Figure 3. From top to bottom: proton kinetic en- 
ergy, acollinearity and missing transverse momen- 
tum for the golden sub-sample. The continuous 
red (green) line is the expectation from FLUKA 
(Saxon- Woods) convoluted with the detector re- 
sponse. 



mation of the asymptotic states through nuclear 
evaporation, hadronic re-scattering or hadronic 
re-absorption. Several kinematic variables are 
only marginally affected by nuclear matters; in 
this case, the corresponding distributions can be 
reproduced once the z^-nucleon interaction is cor- 
rected for Fermi motion and Pauli blocking using 
e.g. a Saxon- Woods (SW) potential for the nu- 
cleus^. Clearly, purely leptonic variables belong 
to this category. Fig. [3 shows the distributions 
of muon energy, muon transverse momentum and 
the II — v angle. Deviations from expectations 
are visible only for very low fi momenta and the 
corresponding large scattering angles. Tp (Fig. 01 
top panel) is, in principle, strongly influenced by 
the presence of nuclear matter, but, in practice, 
the shape of the distribution is determined by 
the Tp > 50 MeV cut and the requirement of 
full containment in the fiducial volume. Figs. El 
and (top panel) provide a useful consistency 
check, and demonstrate that MC reproduces the 
kinematic selection performed during the scan- 
ning and analysis of the golden sub-sample. 

Other variables embedding the reconstructed 
kinematics of the protons are sensitive to gen- 
uine nuclear effects and should depart from the 
naive Saxon- Woods description. In particular, we 
analyzed both the acollinearity and the missing 
transverse momentum of the event (Fig. [3 mid- 
dle and bottom panels). The former is defined 
as 



A 



acos 



Pxp Px^ 1 Pyp Pyi^ 
\ {pIp+pIp)(p 



(1) 



Pxp and Pyp being the transverse momentum com- 
ponents of the proton and Pxfj, and the corre- 
sponding quantities for the muon. For pure QE 
scattering on a nucleon the muon and the proton 
has to be back-to-back in the transverse plane so 
that the acollinearity is zero. 

In spite of the limited statistics and the ~ 10% 
contamination, the imaging capability of the TPC 



•^In this model, nucleons are treated as a Fermi degenerate 
gas located in a potential well of the form E(r) = Eq /(I + 
e(' where a is 0.6 fm; Ea = 46 MeV is the maximum 

well depth and R is the nuclear radius (iJ = 3.6 fm for 
Ar) ini. 



6 



allows to establish empirically the inadequacy of 
the SW description of the nucleus. An improved 
determination of the level of agreement between 
the theoretical expectations and the data can be 
obtained from a background subtraction proce- 
dure; this is based on the test samples mentioned 
in Sec 14. II If the QE analysis is applied to the 
test sample where a tt*^ is clearly identified, and 
the presence of the 7r° is then ignored, an arti- 
ficial acoUinearity excess is generated; in partic- 
ular, 36% of the sample has A > 60°. Similar 
results are obtained for p™^^^ where 72% of the 
identified resonance events shows a p^^^^ greater 
than 400 MeV. The corresponding background 
subtracted are shown in Fig^The Kolmogorov 
probabilities for acoUinearity (transverse momen- 
tum) are 0.52 (0.30) for Fluka and 0.003 (0.027) 
for SW. 

5. Conclusions 

We discussed the first exposure of a liquid Ar 
TPC to a multi-GeV neutrino beam. The data 
provided relevant information to experimentally 
establish the effectiveness of the LAr technology 
in the reconstruction of low-multiplicity neutrino 
interactions. In spite of the limited size of the de- 
tector, nuclear effects beyond Fermi motion and 
Pauli blocking have been observed as perturba- 
tions of the quasi-elastic v^, CC interaction kine- 
matics. 
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Figure 4. Background subtracted cumulative dis- 
tribution of missing transverse momentum (top) 
and acoUinearity (bottom). 
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